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Abstract 

We consider the imitative monomer-dimer model on the complete graph introduced in [T]. It 
was understood that this model is described by the monomer density and has a phase transition 
along certain critical line. By reverting the model to a weighted Curie-Weiss model with hard core 
interaction, we establish the complete description of the fluctuation properties of the monomer 
density on the full parameter space via Stein’s method of exchangeable pairs. We show that this 
quantity exhibits the central limit theorem away from the critical line and enjoys a non-normal 
limit theorem at criticality with normalized exponent 3/4. Furthermore, our approach also allows 
to obtain the conditional central limit theorems along the critical line. In all these results, the 
Berry-Esseen inequalities for the Kolomogorov-Smirnov distance are given. 
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1 Introduction and main results 

In [H [2], the authors introduced a mean-field system of interacting monomers and dimers with 
imitative interaction, called the imitative monomer-dimer (IMD) model. Depending on the attractive 
potential J > 0 and the monomer potential h, this model was described by the infinite volume 
limit of the monomer density and was shown to exhibit a phase transition along a critical line 
r. More precisely, the monomer density converges to a constant mo{J,h) for any {J,h) ^ T and is 
concentrated at two distinct values mi{J, h) and m 2 {J, h) for (J, h) G F. The aim of this investigation 
is to establish limit theorems for the monomer density on the full parameter space (J, h). Following 
the methodology of the Curie-Weiss (CW) model [6l[7], the previous known results in this direction 
were obtained recently in [2j, where it was proven that this quantity satisfies the central limit 
theorem away from the critical line and possesses a non-normal limit behavior at the criticality 
with normalized exponent 3/4. In the present paper, we study the limit theorems for the monomer 
density in a completely different approach via Stein’s method for exchangeable pairs. We recover the 
results in [2] and provide the Berry-Esseen type inequalities for the Kolmogorov-Smirnov distance. 
Furthermore, our approach also extends to the conditional central limit theorems along the critical 
line F given the monomer density being above or below any fixed level ^ between m\{J,h) and 
m 2 (J, h). These results together conclude a complete description of the fluctuation properties of the 
IMD model. 

We now introduce the IMD model and state our main results as follows. For V > 1, let C = {V, E) 
be a complete graph with vertex set V = {1,..., N} and edge set E = {uv = {u, v} u,v u < 
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v}. A dimer configuration D on C is a set of edges such that uw ^ D for all w ^ v if uv & D and the 
set of monomers ^{D), associated to -D, is the collection of dimer-free vertices. Denote by S> the 
set of all dimer configurations. Apparently, by definition, the dimer configuration and the monomer 
set satisfy the equation of hard core interaction, 

2\D\ + \^{D)\= N. (1) 

The Hamiltonian of the IMD model with imitation coefficient J > 0 and external field h £ M is 
defined as 


for all D G where 


-H{D) = N{am{Df hm{D)) 

\^m 


m{D) = 


N 


is called the monomer density and the parameters a and b are given by 

T Ah Ai 

a = J and b = — - - \- fi — J. 

The associated Gibbs measure and free energy are defined respectively as 

e-mo) 


F{D) = 




and 


PN = 






It is well-known that the infinite volume limit of the free energy of the IMD model is given by 


where letting 


and 


lim ptv = sup p{m), 
m€[0,l] 


1 


g{x) = - (ve^^H-2e^ — 


( 2 ) 


( 3 ) 


( 4 ) 


r(x) = (2x — 1) J -|- h, 

the function p is defined as 

p{m) = —Jm^ — —(^1 — g o T{m) + log(l — go T{m 

In [1], it has been investigated that the IMD model exhibits three different phases. We summarize 
the first two as follows. Let 

Jc = - i=r- and he = - log(2\/2 — 2) —-. 

4(3-2^2) " 2 ^4 
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There exists a function 7 : (Jc,oo) —>■ M with 7 (Jc) = such that for T := {(J, 7 (J)) : J > Jc}) if 
(J, /i) ^ r, then (l 2 |) has a unique maximizer mo and this quantity satisfies 

mQ= gorimo). (5) 

Furthermore, if (J, h) / (Jc, /ic), then p"(mo) < 0 and if (J, /i) = (Jc, he), then mo = me := 2 — \/2 
and 


p'{mc) = 0, p''{mc) = 0, p^^\mc) = 0, pi^i(mc) < 0. 


The importance of the maximizer lies on the fact that the monomer density satisfies the law of large 
numbers that m{D) mo for any (J, h) ^ T, which can be seen either from [21 Theorem 1.5] or from 
Lemma [2] below. It is therefore natural to investigate the fluctuation of the monomer density, for 
which results related to this direction have been implemented in a recent paper [ 2 ], where the authors 
proved the limit theorems for any pair (J, /i) ^ T by adapting the classical treatment for the CW 
model from PEI- Our main results here establish the same limit theorems and more importantly, 
give the Berry-Esseen type inequalities. 


Theorem 1. If (J, h) ^ T U {(Jc, he)}, then there exists some constant K such that 


sup 

Z 


r \^{D)\-Nmo 
V iVi/2 


<z) - ¥{X < z] 


< 


K 


( 6 ) 


where X is a normal random variable with mean zero and variance A := —p''{mo) ^ — (2J) ^ > 0. 
If {J,h) = {Je,he), then there exists some constant K such that 


sup 


/|.^(T>)| - A^mo 




iV3/4 


< z) 


< 


K 

iVl/4’ 


(7) 


where letting Ac := —p^^\mc) > 0, the random variable Y has density ce Yz‘^I 2 A ^ ^ normalizing 
constant. 


The third phase is along the critical line {J,h) G T. The variational formula (I2|) now has two 
distinct maximizers mi and m 2 with mi < m 2 . They both satisfy ([5]) and p”{mi),p"{m 2 ) < 0, 
see [1]. Moreover, it was later understood in [2] that the monomer density converges weakly to an 
atomic measure pidm^ +P 2 <^m 2 for some pi,P 2 > 0. Let ^ G {mi, m 2 ). Set 

^i = {Dg^: m{D) < 

&2 = {De^: m{D) >i}. ^ ^ 


Our next result presents central limit theorems for the monomer density conditioning on Qi and & 2 - 
Theorem 2. If (J, h) G T, then there exists some K > 0 such that 


sup 


j \.Y^{D)\-Nmi 
V iVi/2 


< z 


- < z) 


< 


K 


where Xg is a normal random variable with mean zero and variance Xi := —p"{me) ^ — ( 2 J) ^ > 0 . 

Theorems [Hand [2] together describe the fluctuation properties of the model on the full parameter 
space. They will be established based on the general framework of the Stein method of exchangeable 
pairs in [3]. This approach has been greatly used in obtaining the limit theorems for the magneti¬ 
zation in the classical CW model PHIS] or the mean-field Heisenberg model p. The idea of our 
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argument is to reformulate the IMD model as a CW model with additional weights described by the 
hard core interaction ©• For A'" > 1, set S = {0,1}'^. Define a Hamiltonian 

—H{a) = N{am{a)‘^ + bm{a)) 


for any a = (cti, ..., un) G S, where 



i=l 


is called the magnetization of the configuration a. In addition, we denote by s/ (u) the set of all 
sites z G F with fij = 1 and by D[a) the total number of admissible dimer configurations D G ^ 
that satisfy ^{D) = Using these notations, we introduce the Gibbs measure, 


P((t) 


D{a) exp(—i/(cr)) 

exp(-H'(r))' 


(9) 


In other words, this defines a weighted CW model on S and more importantly, from the following 
identity 


^ 1(|^((t)| = t)D{a)exp{-H{a)) = ^1(|^(D)| = t)exp{-H{D)) 

a D 

for any t = 0,1 ,... , A^, it satisfies 

P(m(a) = l) =p(m(D) = l). (10) 

From this equation, to prove the limit theorems for the monomer density in the IMD model, it suffices 
to investigate the magnetization in the weighted CW model. We remark that while the space of 
all admissible dimer configurations & is not a product space, the hypercube S has a nice product 
structure. The main difficulty one needs to handle here is the effects brought by the weights D{a). 
Indeed, following the equation of the hard core interaction m, they are equal to zero if Af — (cj)| 

is not even and they are given by a large combinatorial number (see (jI9p ) if A^ — |j 2 /(cr)| is even. 
Recall the scheme of Stein’s method to establishing the limit theorems for the magnetization in the 
classical CW model 13111115]. One constructs the exchangeable pair for the sampled configuration a 
by choosing a site i uniformly at random from V and then replacing ai by (t(, whose law follows the 
conditional distribution of a given {aj)j^i and is independent of Uj. The key step is the conditional 
moment computations for the exchangeable pair, where the resulting equations are clearly related 
to the derivatives of the infinite volume limit of the free energy. In the present case, due to weights 
D{a), we shall construct the exchangeable pair for the sampled configuration a from the Gibbs 
measure (111) by updating a pair of spins {ai,aj) at a time rather than just a single spin. As one shall 
see, the relations between the corresponding conditional moment computations of the exchangeable 
pair and the derivatives of the infinite volume limit of the free energy now become very indirect. 
Several derivations are non-typical and more subtle computations are needed compared to those in 
the CW model. 

Acknowledgements. The author thanks Pierluigi Contucci for several enlightening discussions on 
the monomer-dimer model and bringing the results in [2| to his attention, which lead to the current 
work. The author is indebted to Qi-Man Shao for the fruitful discussions about the Stein’s method 
and to the Department of Statistics in the Chinese University of Hong Kong for the hospitality 
during his visit. This research is supported by NSF grant DMS-1513605 and NSF-Simon Travel 
Grant 2014-2015. 
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2 Stein’s method 

In this section, we describe the formulation of Stein’s method from [3]. The exchangeable pairs we 
shall use in this paper will be constructed through the following general proposition. 

Proposition 1. Let C = iCij ■ ■ ■ iCn) be a N-dimensional random veetor and uv be sampled uni¬ 
formly at random from E. Let {CuXv) be the conditional distribution of {(uXv) given {Ci)i^u,v and 
be independent of {Cu,Cv)- Set (' = ) • • • ) Civ) with (' = (i for all i u,v. Then ( and (' are 
exchangeable. 

Proof Let F, G be any bounded measurable functions and denote by the cr-algebra gener¬ 
ated by (Ci)i^u,v Note that conditioning on ^uv, ((u,Cv) and (Cv,C^) are i.i.d. Consequently, the 
exchangeablility of f and C' follows by 

uv£E 

= ^ E E[E[FiC)\,^uv]E[G{C'Wuv]] 

uvGE 

= ^ E E[E[Fia\.^uv]E[GiCW^,]] 

uv£E 

= E[F{C')G{C)]. 


□ 

Suppose that g{t) is nondecreasing and g{t) > 0 for t > 0 and g(t) < 0 for t < 0. Let Z be a 
random variable with density 

p{t) = cie-^ ^0 

for t € R, where cq > 0 and ci is the normalizing constant. Let A = W — W'. Then Stein’s method 
for exchange pair yields the following Berry-Esseen type inequality. 

Theorem 3 (Theorem 1.2 [3]). Let (IT, IT') be an exchangeable pair. Assume that there exist two 
real-valued functions g and r on R such that 

E[W-W'\W]= g{W) + r{W), (11) 

Suppose that there exists C 2 < oo such that 

coW{x)\(\x\ + —) min(—, -— yE\) - 
V Cl/ Vci \cog[x)\t 


If \A\ < 6, then 


sup |P(IT <z)- P{Z <z)\<2,El- ^E[A^\W] 


+ —E|r(IT)| 
Cl 


-I-Cl max(l,C2)(I + (5^Co I (^2 + y 


+ 


ClC 2 


(13) 
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3 Moment computations away from the critical line F 

Throughout the rest of the paper, we shall use K to stand for a positive constant that is independent 
of N and could be different at each occurrence. For any t £ T,, uv £ E and s,t = 0, 1, we use the 
notation to denote the configuration /? G S that satisfies pi = Tj for all i ^ u,v and Pu = s and 
Py = t. Let us sample a from P and let uv be sampled uniformly at random from E. We define 
{ay,ay) as the conditional distribution of ((Tu,o'y) given {cri)i^u,v and independent of (au,o'y). In 
other words, 


(o'n = s,a'y = t\a) = 


P(a. 


St \ 
uv) 


P(uii) + P(aiS) + P(uOi) + P(aOO)' 


Note that any dimer configuration D with ^(D) = s^{a) satisfies the equation of the hard core 
interaction, 


N 


2\D\+Y,^, = N, 


2=1 


which deduces that 


D{cr]^y) = D{a^y) =0, if cr„ = = 1 or cry = (Jy = 0, 

D{(y^uy) = D{all) = 0, if au = l,ay = 0 or ay = 0,ay = 1. 


Consequently, if ay = ay = 1, 

= (^v = Ik) = 
= 0-^ = Ok) = 

if 0, 

= Ik) 

= 0 ,o-(, = 0 |o-) 
and if ay = l^ay = 0, 


_ iM _ 

P(o-ii) + P(aOO) D{a) + ’ 

P(a“) _ 2)((j00)g-4am(<7)+4a/iV-2fe 

P(cjii) + Pk2S) ~ D{a) + D(c 700 )e- 4 “™('^)+ 4 a/Af- 2 b ’ 


P((j),k _ D^„^^y^a^(V+ia/N+2b 

P(aik + P(ctOO) “ D{a) + ’ 

_^k)_ 

¥{all) + P(fTOg) D{a) + ’ 


(14) 


(15) 



= Ok) 



1 

P(aiS)+P«) 

D{afy)+D{a^y) 

“ 2 

b 

o' 

II 

= Ik) 


D{a^) 

1 

P(aiS)+P(a0i) 

D{a]^y)+D{a^y) 

2 


(16) 


To sum up, if u, v are either both monomers or both vertices of some dimers, then they only could be 
updated as either both monomers or both vertices of some dimers. If only one of u,v is a monomer, 
then after the update they again has only one monomer. Let a' be the random vector obtained by 
replacing ay and ay by a'y and a'y, respectively. From Proposition [H a and a' are exchangeable. The 
proposition below will play an essential role to control the Berry-Esseen bounds (|13l) . 
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(17) 

(18) 


Proposition 2. Let M = ~ X)£i have that 

E[M — M'\a] = Li(m((T)) + Ri{m{a)), 
E[(M - M')^|cr] = L 2 {m{a)) + R 2 {m{a)), 


where recalling T{m) from 0, Li,L 2 ,Ri, R 2 satisfy that for all m G [0,1], 


Li(m) 

L2{m) 


2(1 — m){m? — (1 — m)e^'^*^™')) 
(1 - m) + e2^(™) 

4(1 — m){mf + (1 — m)e^'^^™^) 
(1 - m) + 


|i?i(m)|, |i? 2 (m)| < ^ 


/or some constant K > 0. 


Proof. Let k = 1,2. Consider 

E[(M - MOV] = E V"k] 

uvGE 

= E E[(o-„ + o-„ - o-v o-^,) V] 

u,v££/[cr):u<v 

+ tE E ^[{(Tu + (T^ - - cr'^f\(T] 

(a) 

+ r|| E m(Tu + (Ty-a'^-al)^\a]. 

U,V^S2/[(7):U<V 


We compute each summation as follows. For u G ^/{a) and v ^ ^(u), since could only be 

either (1,0) or (0,1) from (fTBI) . it follows that 

IE[(u„ + a, - < - a',f\a] = E[(l - l)0u] = 0. 


To compute the first and third summations, note that the total number of L/2 dimers on a complete 
graph of size L can be computed as 


1 (L\ (L -2\ (L- 2(L/2 - 1)\ ^ L\ 

{L/2)l\2)\ 2 J \ 2 ) (L/2)! 


(19) 


Take L = N — \£^{a)\. For u,v G £/(a) with u < v, since (u^,,cr^,) could only be either (1,1) or (0,0) 
from (I14p . this yields 

IE[(fT„ + cr^ - ctV fT/)^|cT] = 2 ^P((t(, = 0,0-/ = 0|(t) 

2fc^(^00)g—4am((T)+4a/Af—2fe 

~ D{a) + L>(cr00)e-4“”^W+4«W-2fe 

r,k (.^^+2)! r)—(L/2+l) „—4:am(a-)+4a/N—2b 
_ ^ (L/2+1)!^ __ 

L\ n-L/2 I L+2)! n-(L/2-^-^)^-4nm.(n^-^-4nJN-2b 

JLJW- + (L/ 2 + 1 )!^ W -r ;e V / 

2^(L + -^'j^—4am,{cr)+4a/N—2b 

]_^g—4am((T)+4a/Af—26 ’ 
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while for u,v ^ £/{a) with u < v, since {a'^,a[,) could only be either (1,1) and (0,0) from (fTCI) . we 
conclude 


E[(ct« +crv-cr'^- (y'vfW] = = 1, cr(, = 1 |cj) 

~ D{a) +D(o-ii)e4“«^W+4a/iV+2fe 

C r)\k (-^~2)! rj—fi^/2—1) ^4amfcr)+4a/A^+26 
_ {L/2-iy, _^_ 

“ ^2-^2 + ^iI0^2-W2-l)e^am{a)+4a/N+2b 
^_2^^g4ccm(f7)+4a/A^H-2b 

_1) H“ * 

Combining these two equations together, we obtain 
E[(M-M')^|u] 

1 //|J^(C7)|\ 2^(L + /|^(o-)'=|\ (^_2)k^4am(a)+4a/N+2b 

l-El I \ 2 y 1 + + l)e-4am(o-)+4a/Ar-26 2 J — I'j Q4:am{a)+Aa/N+2b 

_ 2^-^ /M{a){M{a) - 1)(L + l)e-4am(a)+4a/iV-2fe (_l)fcL(L - l)e4am(-T)+4a/7V+26 
|J5;| V 1 _|_ _|_ X)g-4am((T)+4a/iV-2fe g4am((T)+4a/iV+2fe j 

_ 2^-iiV2 m(cj)(m(cj) - l/iV)(l - m{a) + i/Ar)g-4am(<T)+4a/iV-2fe 

“ |£;| 1/iV + (1 - m((T) + l/iV)e-4“M^)+4a/Af-2b 

(_l)fc2*^-lAr (1 - m(a))(l - m(cT) - l/Ar)g4amH+4a/7V+2b 

''' |E| (1 - m(cr) - 1/A^) + e4“™W+4a/Af+26/jV 


Substituting a = J, _ 7 Vg 2 ft. 2 J — l)/2 into this equation gives 


E[(M-M')*^|o-] 


2^ m{a){m{a) - 1/N){1 - m{a) + /N 

(1 - 1/iV) 1/N + (1 - m{a) + l/Ar)e-2^(™(^))+4J/^/Ar 

(-2)*^ (1 - m(u))(l - m(cT) - 1/Ar)g2r(m(<7))+4J/7V^ 

^ iV(l - 1/iV) (1 - m{a) -1/N) + e2r(m(-T))+4j/iV 

Uk{m{a), 1/iV), 


where for 0 < m < 1 and small t, 


Note that 


2^ m{m — t){l — m + t)e 2T(m)+4Jt 
k{m, C 1 ('l _ ^ i)g-2r(m)+4Jt 

(—2)^ (1 — m)(l — m — 

(1 - t) (1 - m - t) + e2'rM+4'^* 


17fc(m,0) 


2fc^2(i _ ^)g-2r(m) ^ (- 2)^{1 - 

1 + (1 - m)e-2^M (1 - m) + 

2 ^(m^(l — m) + (—1)*^(1 — 

(1 - m) + e2^(™) 

Lk{m). 
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Letting Rk{m) = dtUk{m,t)dt, we have by the fundamental theorem of calculus, 

E[(M - M')’^\a] = Lk{m{a)) + Rk{m{a)), 

where since the numerators in Uk stays away from zero, there exists some K such that |i?fc(m)| < 
K/N. This finishes our proof. 

□ 


4 Proof of Theorem [T] 


Suppose that mo is the unique maximizer of ([2]). Recall M,M' from Proposition [2l For A: = 0 or 1, 
we set 

TT/ _ ^ ~ ^'^0 J TT^/ _ ~ ^^0 

~ JY(2k+l)/(2k+2) “ JY(2k+l)/(2k+2) ■ 

From the previous section, it is easy to see that (W, W') is exchangeable. The following lemma is 
the central ingredient of our argument. 

Lemma 1. Suppose that mo is the unique maximizer of m and that for some integer k = 0 or I, 

Li\mo) = 0 

for all 0 < £ < 2k and 

> 0 . 


We have that 


E[1T - W'\W] = g{W) + r{W), 


where 


r(2fc+l)/ N 

g{W) = ^3.i + ^ 


{2k + l)!iV{2fc+i)/(fc+i) 


and r is the remainder term satisfying 

\r{W)\ < 


K 


j\r(4fc+3)/(2fc+2) 


+ 1 


In addition, 


( 20 ) 


( 21 ) 


( 22 ) 


where 


E 


l-^W.[{W -W'f\W] 


< K 


1 ^ 1 
iVl/(2fc+2) + Jf 


2]^{2k+l)/(k+l) 

L 2 (mo) 


(23) 


(24) 


It should be pointed out that in the CW model mm the exchangeable pairs were constructed by 
choosing a single site i uniformly at random and updating Uj by a{, whose law follows the conditional 
distribution of CTj given {aj)j^i. In those cases, the function g can be expressed in terms of the second 
or fourth derivative of the infinite volume limit of the free energy, but this is not the case in the 
IMD model. 
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Proof of LemmaUi From the given assumptions, the Taylor formula yields 


r I I n I i ,2t+i I XiT WWtr) 

( 2 ^, + !)! +-( 2 FTT)!-■ 


Since 


IT , W' 

m(a) -mo = jyj^/( 2 k+ 2 ) ) - "iQ = 


JYl/(2k+2) ’ 


we have that from (jl7h . 


- ‘’"in'! = = 9(tr) + r(ir), 


where g is given by (I2ip and 


r{W) = 




{2k + l)!iV(2fc+l)/(2A:+2) 


+ 


1 


Ri 


w 


Here 


follows by 


|r(lT)| < K 


(m(cr) — 


+ 


]Y(2k+l)/(2k+2) V Arl/(2fc+2) 


K 


+ mo . 


.(2A; + l)!iV(2fc+l)/(2fc+2) JY(2k+l)/(2k+2) JSfJ JV(4fc+3)/(2fc+2) 

To show (1231) . we use (|18l) and the fundamental theorem of calculus to obtain 


W 


2k+2 


niw - w'?\w] = 

vv ) \VY \ j^^2k+l)/{k+l) ^ ]\j{2k+l)/{k+l) 


9 rm{a) 

- ^ 

Co J mo 


L'2{s) 


and therefore, 
E 


1 - jE[(W - W'f\W] 


= —E 


j\j'{2k+l)/(k-\-l) 


ds “h 


R2{m{a)) 

J\[(2k+l)/{k+l) 


2 

< K 

< K 


rm{a) 


L'2{s) 


'mo 


^^( 2 ^+ 1 )/ (fc+ 1 ) 


ds + 


R2{m{a)) 


J\^{2k-\-l)/ (fc+l) 


(E\m{a) -mo\ + — 
E|TT| , 1 

Ari/(2fc+2) + JfJ- 


+ 1 • 


□ 

Next we need an auxiliary lemma. Denote by O the collection of all maximizers to ([2]). Note 
that from the introduction O contains at most two points. For m E [0,1], let d{m) be the distance 
from m to the set O. 

Lemma 2. For any <5 > 0, there exists g > 0 such that 

F{d{m{a)) > 6) < iFe-^T 
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Proof. Let U = {m G [0,1] : d{m) > 5}. By the virtue of (fTII|l . it suffices to prove that for any 5 > 0, 
there exists r] > 0 such that 

F{m{D) eU)< Ke-^^. 

Note that 

^logP(m(D) G [/) < ^log exp(-i/(D)) - ^log^exp(-iL(D)). 

D:m.{D)eU D 

Set A = {0, l/N ,..., {N — Pj/N, 1}. Observe that 

<^m(D),mexp(-iL(L>)) = 6m{D),m^^vN{am{Df + hm{D)) 

= 6m{D),m exp N{a{2m{D)m — wf ) + bm{D)) 

= S^(^£)j„^expN{{2am + h)m{D) — amf). 

We obtain 


Y exp(-L/'(D)) = ^l(m(L>) G ?7) exp(-iL(D)) 

D:m{D)eU D 

= Y '^^rn{D),mew{-H{D)) 
meAnU D 

= E E exp N{{2am + b)m{D) — amf) 

meAnU D 

< (Ai + 1) sup Y] exp N{2am + b)m{D) 

niGAnU 


and thus, 

1 


— logP((m((T) G U) 


< 


Here, 


log(jV + 1) g^p|_Q ,^2 + ^iog^exp(iV(2am + h)m{D))^ - ^ log^exp(-iL(L>)). 

m€U ^ ^ 

^ log exp {N{2am + b)m{D)) 


is indeed the free energy of an IMD model with parameter (J', h') = (0, (2m — 1) J + h) and its 
thermodynamic limit, according to the formula (I5|), is equal to 

' 1 — o T{m) 


+ log(l — go T{m 


As a consequence. 


limsup — logP(m(Ll) G C/) < sup p{m) — sup p{m) =: —2rj. 

Af—>00 A* m£U me[ 0 ,l] 

Since p{m) is continuous on [0,1] and any point in O is away from the maximizers with distance at 
least 6, we conclude that p > 0 and that there exists some Nq such that for all N > Nq, 

■^log¥{m{D) G [/) < —ry 


and consequently, F{m{D) G H) < Ke 
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Lemma 3. Under the assumptions of Proposition\^ there exists some K > 0 such that < K 

for all N > 1. 


Proof. From (I20p . we have that 


Multiplying W on both sides and then taking expectation give 


EkF 


2fc+2 


(2A: + l)!iV( + )/( + 


L 


( 2 fc+l) , 


1 (”^o) 


< 


{2k + l)\ 


L 


(2k+l) 

1 


{mo) 




iFE|Ty|2^+3 iFE|PF| \ 

jVi/{2fc+2) + j\ri/{2fc+2 )) ’ (25) 


where we have used (I22p to bound r(ld^). Here from the exchangeablility of W and W, we can 
express E[(VF — W')W] = 2“^E(VF — VF')^, which combined with the trivial bound |VF — VF'| < 
2 /A^( 2 ^+i)/( 2 fe+ 2 ) allows to control the first term of (l25p . 

^(2fc+i)/(fe+i)jgj^^ _ ^r')w] < 2. 

As for the third term, we use the bond |VF| < A'i/(2^+2) obtain < 1. To bound 

the second term, for any <5 > 0, Lemma [2] says that there exists some p > 0 and K > 0 such that 

P(|VF| > J7VV(2fc+2)) ^ P(|m(f 7 ) - mol > < 5 ) < Ke-^^ 


for all > 1 . Consequently, using again the trivial bound |VF| < A^i/( 2 *:+ 2 )^ 

E|TT|2fc+3 _ E[|fF|2^+3; |fF| < 5iVl/(2A:+2)] ^ E[|kF|2*'+3; |kF| > 5iVV(2fc+2)] 

iVV(2fc+2) “ j\rl/(2fc+2) ^ ^l/(2fc+2) 

< 6 E\W\‘^^+^ + lVP(|m(a) - mo| > 5 ) 

< 6E\W\^^+^ + KNe-^^. 


Plugging these three bounds into (l25P gives 


1 - 


A:( 2 fc + 1 )! 

Lf+')(mo)' 


2fe+2 


< 


(2fc + l)! 

Lf+')(mo) 


(^2 + KNe-^^ + , 


which completes our proof. 


(26) 


Lemma 4. Suppose that the conditions of Proposition\^hold. Let Z be a continuous random variable 
on M with density 


p{z) = Cl exp(^-d 2 ;^^+^^. 


where 




( 2 fc+l) 


{mo) 


{2k + 2 )!L 2 (mo) ’ 

and Cl is a normalizing constant. Then there exists some constant K independent of N such that 

K 


sup 


P(kF < z)-F{Z < z) 


< 


JY 1 /( 2 A:+ 2 ) ' 
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Proof. Recall cq from (|24p and g,r from (1201) . We define 

pit) = cie-^^'^3{s)ds ^ 

where ci is a normalizing constant such that p is a probability density on M. Using these cq, ci,g, r, 
we now verify m for some C 2 , which can be easily seen since 

= (2k + (|x| + 4) mm(i 

has a limit at infinity and is clearly bounded for arbitrary small x. As a result, the inequality (jl3p 
with (5 = 2 /Ar( 2 fc+i)/( 2 fc+ 2 ) gives 


sup |P(VU < z) - F(Z < z)\ 


< 3K 


1 


+ 


< 


]yl/(2k+2) 

2 ci max(l, C 2 ) 
jY(2A:+l)/(2fc+2) 

K 




+ 


AK 


]Sfl/{2k+2) 
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(E|lUp ^+2 ^ 

C 2 d„ 


_/\r(2fc+l)/(2A:+2)^2(mo) 


{(2 + ^E|lU| 


2fc+l 


+ 


C 1 C 2 


jSll/{2k+2) ’ 

where the first inequality used Lemma [1] and the second one used Lemma [3l 


□ 

Proof of TheoremUi Recall A, Ac from Theorem [H Suppose that (J, h) ^ T U {(Jc,/ic)} and mo is 
the unique maximizer of (l2|). From ([3|), mo satisfies 


2 mo + 


(27) 


or equivalently 

mo = (1 — mo)e^'^*'™'°^ 


(28) 


Note that from ([^ . 


p"(mo) = 2J{2Jg’ o r(mo) - 1) 

g4T(mo) _j_ 2g2r(mo) 

— 2J\ 2J 


_ g2T(mo)\ 


' ^g4T(mo) _|_ 4g2T(mo) 

Oj2mo + (4J(mo - 1) + 

~ 2mn + e2^(™o) 


- 1 


and thus, 


A 



1 

p"(mo) 



2(1 - mo)e2^(”"o) 

2mo + (4J(mo - 1) + l)e‘^prao )' 
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On the other hand, (f25|l implies that Li(mo) = 0 and that 


2L[{mo) 1 2mo + (4J(mo — 1) + 2mo + (4J(mo — 1) + 1 

2!L2(mo) 2 mg + (1 — 7710 ) 6 ^"^ 4(1 — mo)e^'^(”*o) 2A 

Since the equation (l?7l) also implies 


p"{mo) + 2J = 2Jg' o T{mo) 


4J(1 - 

2mo + e2^(”^o) 


we conclude that A > 0 and thus, L[{mQ) > 0. Lemma [Hand (fTOl) then deduce Q. Next assume 
that (J, h) = {Jc, he). In this case mo = rric = 2 — y/2 and a direct computation gives 


Li{mc) = 0 , L[{mc) = 0 , L'{{mc) = 0 , L'”{mc) = 6 + 


17^2 


and 

2Lf (me) _ 1 17\/2 _ p^^)(mo) _ ^ 
4!L2(mc) “ 2 48 “ 24 “24' 

Lemma [Hand (1101) then yield (|7]). This finishes our proof. 


□ 


5 Moment computations along the critical line F 


This section is devoted to dealing with some moment computations for the parameters along the 
critical line. Set M(p) = N~^ Pi p = (pi,, pjsi) G S. Let 

Si = [pe^-.M{p)<iN], 

^2 = {/5 G S : M(/5) > ^iv}. 

For i>l, define the probability measure = P( • on S^. For the same reason as (fTOjl . one sees 
that 



where is defined in ([ 8 |). Therefore, to prove the conditional central limit theorem for the monomer 
density in Theorem |2l it suffices to establish the central limit theorem for the magnetization under 
P^. Following the same construction as before, let a be sampled from P^ and uv be a uniform random 
variable from E. Under the probability measure P 7 , let (cjF< 7 (,) be the conditional distribution of 
(cr^, fj.,;) given {ai)i^u,v and independent of (cj.i;,cj^), that is. 


¥i{a'u = s,cr(, = tier) 





This pair {cF,a') is therefore exchangeable from Proposition [TJ Set 


= {u : M(cj) < CN- 2 Y 
5" = |u : - 2 < M{a) < 
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and 


5^ = |cj : M(cj) > eA^ + 2|, 

5" = |cj : ^iV + 2 > M(cj) > ^Ny 


In the case a G S'^, we have 
= s,a'^ = t\a) = 


n<.)/nsi) 


= F{a'^ = s,a[, = t\a), (30) 


(P(aii) +P(aiO) + P(a0i) +P(a00))/P(5,) 
where we used <7^1,, ^ However, if ct E S", this equation is no longer valid. 

Lemma 5. Recall Lk and Rk from Proposition\^ We have that 

Ei[{M - M')^|M] = Lk{m{a)) + Rk{m{cj)) + r^,fc(m(a)), 
where 'K(\T(,^k{fn{a))\ < K¥ii{a E S”). 

Proof. Observe that if cr E S"^, from (l30]l and then (flTl) . (fTHI) . 

Ee[{M - M')’^\a] = ^ E 

uvdE 


(31) 


Thus, 


uvGE 

= K[{M - M'f\a] 

= Lk{m{a)) + Rk{m{a)). 

E(,[{M - M'f\a] = l{a E S'^){Lk{m{a)) + Rk{m{a))) + l{a E S';)Ee[{M - M')’^\a] 
= Lk{m{a)) + Rkim{a)) + T{m{a)), 


where 

T{m{cr)) = 1((T E S") {Ei[{M - M')^\a] - Lk{m{a)) - Rk{m{a))). 

Taking conditional expectation E£[-|M] and letting Ti^k{m{a)) = E£[T{m{a))\M] give ([31]) . Here 
Ei\T^^k{i^i(^))\ < K^eio' E S”) holds true by applying the trivial bound \M — M'\ < 2 and the fact 
that Lk,Rk are bounded. 

□ 


6 Proof of Theorem [2] 


In this section, we suppose that (J, h) E T. This assumption implies that p'{m£) = 0 and p''{m() < 0 
and therefore using the first half of the derivation of the proof of Theorem [H they can be transferred 
as Li{m(:) = 0 and L[{m£) > 0. Furthermore, the quantity defined in Theorem [2] is positive. 
Denote 


W^ = 


M — Nrrii 
A^V2 


and = 


M' — N mi 
]VV 2 


First by adapting exactly the same argument as the proof of Lemma [T] and applying Lemma [5l one 
obtains an analogue of Lemma [T] for the conditional probability P^, whose proof will be omitted. 
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Lemma 6. We have 


MW I - w',\Wi\ = g^iw,) + n{w^), 


(32) 


where 


and r£ satisfies 


In addition, 


St(We) = 




Ef 


CO,£t 


i-^M{Wi-wM\Wi\ 


< K 


1 


+ + iMM ”)), 


iVi/2 ' N N 


where cq^i = 2 N/L 2 {mfi). 

The next lemma below will play a similar role as Lemma [3l 
Lemma 7. We have that E^lW^p < K for all N > 1. 

Proof. Similar to the proof of Lemma [3l multiplying Wi on both sides of (I32h and then taking 
expectation give 

^iWl = (EeiiWe - Wl,)We] - E^W^We)) 

(iVE,[(TT, - W^)We] + G S'f)). (33) 


< 


L' 


Let us now bound each term on the right-hand side as follows. Using the exchangeability of (W£, W^) 
under and the bound SfWn — fU^'l < 2/A^^/^, we obtain the control of the first term 


NEfiW^ - W'fifWi = yE^(lUf - W'M < 2. 


For the third term, it can be easily controlled 

Ari /2 


< 1 


(34) 


(35) 


by noting IVF^I < N^/'^. As for the other two terms, note that since P((T G S^) converges to > 0 as 
N ^ oo, this together with Lemma [2] implies that for any 5 > 0, there exists some g > 0 and A' > 0 
such that 

P.dmW -m.l > i) < - “.I a < a < 

r[bi) 


P2(|m((T) — m2| > 5) < 


\m{a) — 7712 ! > (5, m{a) > ff) 

PGS^) 


< Ke 


— 7]N 


for all > 1. Consequently, these imply that for N large enough, 

Fi{a G S'l) < P£(|m(cj) - ^1 < 2/N) < 


(36) 
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(37) 


Finally, proceeding in the same way as (ESI), we see that 

+ KNe-^^. 

Plugging dMI), ([35]), ([36|) and ([37]) into ([33]) gives 

(' - s lifs) ('+’ 

for some K' > 0. Note that K is independent of 5. We can choose 5 small enough so that the 
announced statement holds. 


Proof of Theorem O Recall co,£ and gi from ([6]) . Define 

p^{t) = ^0 


where ci^g is the normalizing constant such that pi forms a probability density on M and di = 
L'i{mi)/L 2 {mi) > 0. Observe that 

co,i\9e{x)\(\x\ + —) min(—, 

V ci,f/ Vci,^ \co,£ge{x)\y 

= 2de\x\(\x\ + —) min( — 

^ Cl,£7 Vci,£ dex^/ 

has a limit at infinity and is clearly bounded for arbitrary small x. This function has a uniform 
upper bound over M that is denoted by C 2 ,£. Let 5^ = 2/A^^/^ and = Wi — Wf Applying these 
co,i,ci/,C 2 /,ge,r£,pi and {W£,W^) under the inequality ifTT]) leads to 


sup|P^(W,<z)-P(A^<z)| 


< 3E, 

+ 

< K 


1 - ^E,[A2|W,] 


+ ^E,|r,(lP£)| 

Cl,£ 


ci,£max(l,C2/)E£|(5£| + l^^l^co/j (^2 + ^E£|co,£5'£(iL£)|) + 

(^E,|»'|| +1 + ^ + ^ + e S'l) + ip,(, e ST)) . 


Here the second inequality used Lemmas [5] and [6] and |Af| < Finally, using Lemma [71 the 
inequality (|36h and the identity (j29h finishes our proof. 

□ 
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